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ABSTRACT 

We extend recent work that included the effect of pressure forces to derive the pre- 
cession rate of eccentric accretion discs in cataclysmic variables to the case of double 
degenerate systems. We find that the logical scaling of the pressure force in such sys- 
tems results in predictions of unrealistically high primary masses. Using the prototype 
AM CVn as a calibrator for the magnitude of the effect, we find that there is no scal- 
ing that applies consistently to all the systems in the class. We discuss the reasons for 
the lack of a superhump period to mass ratio relationship analogous to that known 
for SU UMa systems and suggest that this is because these secondaries do not have 
a single valued mass-radius relationship. We highlight the unreliability of mass-ratios 
derived by applying the SU UMa expression to the AM CVn binaries. 

Key words: stars: binaries: close - stars: novae: cataclysmic variables - accretion, 
accretion discs 



1 INTRODUCTION 

The standard model of a cataclysmic variable (CV) is of a semi-detached binary consisting of a white dwarf accreting material 
through an accretion disc from a late-type main-sequence secondary. Among the taxonomical classes of CV are the dwarf 
novae (DNe) and their sub-class named after the prototype SU UMa. DNe show semi-regular outbursts of increased luminosity 
(typically 2-5 mag.) lasting for a few days and recurring on timescales of tens of days to tens of years. The SU UMa systems 
show occasional superoutbursts with a brighter maximum (~ 0.7 mag.) and longer duration (approximately 5 times) than 
normal outbursts. They also show "superhumps" : a periodicity o bserved in the lightcurve during an outburst occurring on a 
period (P s h) a few percent longer than the orbital period (P or b) (|Warnerlll995h . 

DN outbursts are believed to result from an ionization instability that sets in when the accretion disc exceeds a critical 
maximum surface density (E max ). On crossing this threshold, the disc switches to a "hot" state with a significantly higher 
viscosity and mass throughput. This results in a decrease in the surface density until the disc reaches a second critical surface 
density (E m i n ) where the disc can no longer mainta in the hot state. At this point, the disc switches ba ck to a low viscosity 



configuration, begins to refill and the cycle repeats (jMever fc Mever-Hofmeisterlll98ll ; iFrank et al.l ll992'). 



AM CVn systems are analogues of DNe where the mass-donating secondary is a helium-rich degenerate or semi-degenerate 
star. These systems can similarly exhibit outbursts with detectable superhumps. 

Superhum ps arise from an el lipticity induced in the shape of the accretion disc by tidal excitation of a resonance with 
the secondary 1 Whitehurst 1988). This elliptical pattern has a precession period (P pr ) much longer than P or b- Thus the 



secondary and the major axis of the disc come into the same relative alignment on the beat period between these two which 
is consequently just slight ly longe r than P o r b- It is the modulation of the tidal dissipation on this beat period that is picked 
up observationally as P sh . iLubowl l|l991al lbL Il992h derived the steady state precession rate ui pr for an eccentric disc as: 

UJpr — Wdyn + impress (1) 

where Udyn is the dynamical precession frequency and a; prcss is a pressure related term. These authors recognised the similarity 
between the precession of the accretion disc and the familiar inner Lindblad resonance that produces spiral waves in galactic 
discs. Historically, only the dynamical frequency 
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( Hirose fc Osakilll990h las been considered important where q = M%/Mi is the mass ratio, r is the distance from the primary 
(expressed here as a fraction of the separation d) and 2F1 is the hypergeometric function. Recently the necessity of including 

(3) 
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has been reasserted (|Murravll20o"ol ). where u p is the orbital angular frequency of disc material at a radius r (hereaft e r mea sured 
in dimensional units) and i is the pitch ang l e of th e induced spiral wave. This was studied bv lGoodchild fc Qgilviel (2OO6) with 
a detailed integration scheme and Pearson (2006) with an algebraic approach to derive an improved relationship between ui pT 



and the system parameters. In particular, the latter paper used the observed values of P s h to derive values for the component 
masses of 88 CVs. These papers provi ded a theoretical foun dation for the observed correlation between the observed period 
excess (e = (P s h — -Porb)/-Porb) and q ( Patterson et al. 20051 '). While the nature of the superhumps in ordinary CVs is well- 
established, the application of the theory to AM CVn systems relies on a "by analogy" argument. Here we examine the degree 
to which this same approach can be successfully applied to AM CVn systems. 



2 PRECESSION FREQUENCY 
2.1 Dynamical Precession 

The dynamical precession rate can be evaluated by idealising the response of the disc as that of a ring of material at a radius 
r = j~3(l + g) _3 d that has an orbital frequency about the primary in a resonance with the secondary. The expression given 
in ((2} can be then be reduced to 



UJdyn _ O q \ ~* On 



where 
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( Pearson! liooj ) and j = 3 is the appropriate resonance. If the dynamical term were the only contribution to the precession 
rate we could rewrite equation @ as 



2.57i^i P orb 



(6) 



by evaluating the summation with a " typica l" q = 0.16 (cp. Warner ( 19951 )). The coefficient in ((6| differs from the often used 
value of « 3.85. As noted by I Murray! 1 20001 ). this latter value results from an erroneous factor of 2/3 being introduced. This 
comes from considering the time for a test-particle to return to the same relative alignment with the secondary rather than 
a ring of material along the resonant orbit. 

The above prescription for the radius assumes that u) pr <C <^ or b in order to reduce the resonance condition 

U - !)( W P - ^pr) = j(Wp - ^orb) (7) 

to the approximate u) v ~ jc^orb- As such, the infinite sum in ((4} does not carry the accuracy that might normally be assumed 
for such a summation. It also tacitly assumes that the resonance can be represented by the response at a single radius rather 
than spread throughout the disc structure. In any case, we also need to account for the additional pressure related term. 



2.2 General Form of uo n 



A fuller explanation of the following derivation is given in lPearson (2006) but we take this opportunity to correct the numerical 
value of constants given there and to write the equations in a form applicable to both hydrogen- and helium-dominated discs. 
For a disc opacity law 

K = K pT~ 3 - 5 (8) 
where 



Ko = 2.8 x 10 20 m 5 kg -2 K -3 ' 5 



(9) 



is appropriate ( Cannizzo 1992b ). equation A3 of Cannizzo fc Reifi 1 19921 ) can be written as 
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where Eo,h = 406 kg m~ 2 , M\ is the primary mass in solar units and M the mass transport rate through the disc. The 
subscripts h ( c ) will be used to refer to the disc in its hot (cold) state and the subscripts h and He will be used for hydrogen- 
and helium-dominated discs respectively. Assuming free- free opacity is dominant, the relationship 

Ko.ff oc (1 + X)(1-Z) (11) 

( Bowers fc Deeminelll984h allows us to scale to S , H c = 428 kg m" 2 by using X n = 0.7, Y H = 0.27, Z H = 0.03, A Hc = 0.0, 
y He = 0.97, Z He = 0.03 jTsugawa fc Osakill 19971 ). 



Cannizzo. Shaffer fc Wheeler 1 19881 ) give the maximum possible surface density in the cold state as 



-2 ( r \ 105 
E max , H = 114 kg m ( Yo^ ) Ml " 



0.35 -0.86 



for a H-dominated disc and similarly iTsugawa fc Osakil fl997) 



E max , Hc = 2510 kg m -2 (y^— ) Mr°' 33 ( 



0.3 



(12) 
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for one that is He-dominated. For consistency in the form of the equations, we scale the value of this latter relationship to 
a c = 1 and adopt here, 

S max = E crit f -£—) 1,05 Mr°- 35 a-°- 86 (14) 



where £ C rit,H = 114 kg m -2 and £ C rit,Hc = 1080 kg m -2 . 

Integrating ()10|) with respect to radius gives the disc mass Md and similarly (| 14f) gives the maximum disc mass in the 
cold state M max . If the disc fills to a constant fraction / of its maximum possible mass, the equality Md = /M max results in 

f \ l-« f n s -1.23 



M = Mi 
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where M , H = 2.62 x 1 0~ s kg s" 1 m" 2 57 and M ,H e = 60.3 x 10" 8 kg s _1 m" 



Cannizzo fc Reifj jl992l ) gives 



where T , H = 0.246 K m~ o m s5 and T , H c = 0.597 K m" ' 771 si. Hence, 
2 ikT 



/i h m H 



/ n, \ 0. 142 / n s -0.369 / f \ l 
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(15) 

Substituting this into equation Al of 

(16) 

(17) 
(18) 



where Cq iH = 3.40 x 10 3 m 1 ' 229 s 2 and c 2 , Hc = 8.27 x 10 3 m 1 229 s 2. Putting this into ©, we can further make further 
eliminations: 



where 

0.49g"§ 



Eiq- 1 ) = - 

0.Qq~z +ln(l + q~i) 

( Eggleton Ill983l ). The resonance requires 



3worb = 3 {GM 1 (1 + q) Af }3 d~ 5 



(19) 



(20) 



(21) 



and 

r = 3 _ * (1 + q)~? d. (22) 
Equating the radius of the secondary R2 to the size of its Roche lobe Rl,2 ~ E(q)d, we can also substitute for the separation 
R2 



d 



(23) 



The final expression for the pressure term is thus 
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,,, . /^s 0.142 / N -0.369 / f \ 0.429 , o N 0.771 p/ -1\0.771 

where 770, h = 1-27 x 10~ 12 m -1 ' 021 and 770, He = 3.08 x 10~ 12 m 021 . To make further progress we need a form for the 
mass-radius relation i?2(M 2 ). 

2.2.1 Main Sequence Secondary 

Using the main sequence mass-radius relationship 

R 2 = 0.91M 2 °- 75 i? Q (25) 



I Smith fc Dhillonlll998l ). (|24)| can be written as 



= ~3 6 — t/h (26) 

Worb £(g) 1021 (l + g)T5 

where 

/ rv,, \ 0.142 / s -0.369 / f \ 0A29 / R \ °- 771 

,„ = i, 3 *io-w,( S ) (,*) (A) „,-*-. m 

Assuming r/n to be a constant (r? H ), Pearson 1 20061 ) found an excellent best fit value of fj H — 0.0109 for those superhumping 



CVs with independently determined values of q. The summation of §4$ and (|24p results in an algebraically unwieldy, but 
numerically easy to invert, expression for ^^-(g). 

While a constant 7/ H produced an excellent fit to the data, an accurate calculation of the effect in individual systems 
should include the Mi dependence of rju ie. 

/ \ -0.242 

m - j Vh (28) 

where Mi = 0.76 is the weighted mean primary mass of the sample used to determine the value of rj H . This implies that the 
coefficient of proportionality between P s h and P or b is no longer strictly a function of q alone (cf. equation (]6])). However, the 
combination of (|21[1 . (|23[) and (|25[) produces a second constraint P or b = Poib(M2,q). We can thus solve simultaneously the 
^orb(M2, q) and P s h(Mi, M2) expressions for both Mi and M2 together. In principle, this will change the value of Mi and the 
expressions should be iterated upon for a fully self-consistent solution. However, for a main- sequence second ary, allowing for 



this dependence produces masses that differ by only a few percent from those in Table 4 of iPearson (2006). The equivalent 
correction for systems with a (semi-)degenerate secondary is more significant. 



2.2.2 Degenerate and Semi-degenerate Secondary 



The approach of the previous section can be applied mutatis mutandis to a secondary with a non-main-sequence structure. 
Equation (|26|l remains applicable but now 



2.67 
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*r- 0.766„ 

M 2 t?h 



0.0127 



Mo 



(29) 



(30) 



using the compositions given earlier to derive /ih.H = 0.618 and /Jh.Ho = 1.347 and the previously found best value of ?7 H . 

Two forms of the mass- radius relation have b een widely used in the literature corresponding to two possible formation 
channels for these systems ( Nelemans et al. 200lh . If the donor is the remnant core of a low- mass helium star that begins 



mass-transfer before helium burning, a zero temperature degenerate structure is appropriate. For a fully degenerate helium 
white dwarf secondary the mass-radius relation is 

R2 



= 0.0106 - 0.0064 In M 2 + 0.0015M 2 



(31) 



1 Zapolskv fc Salpeter|[l969l ; Rappaport &: Jossl 1984 ). Alternatively, the donor star may begin helium burning before the onset 
of mass-transfer. In this case (a "semi-degenerate" secondary), the mass-radius relationship has been approximated by 

R2 



R, 



= bM~ 



(32) 



wherelTutukov fc Fedoroval(|l989h have b — 0.043, a = 0.062 (hereafter TF parameters) and lSavoniie. de Kool fc van den Heuvell 
1 198fj) found b = 0.029, a = 0.19 (hereafter SKH parameters). 
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Table 1. Component masses derived from th e application of the re lationship established for SU UMa systems to all the AM CVn 
systems with measured P or b and P s h listed in iRitter fc Kolbl (|2003h . 



System 


TF (dyn. 


only) 




TF 




Degenerat 


e 




Eggleton 




-forb 


e 




9 


Mi 


M 2 


'/ 


Mi 


M 2 


'/ 


Mi 


M 2 


'/ 


Afi 


M 2 


W 






AM CVn 


0.056 


2.02 


0.113 


0.064 


1.76 


0.114 


0.078 


0.42 


0.033 


0.076 


0.48 


0.036 


1028.8 


1 2 


0.0218 1 


HP Lib 


0.036 


2.72 


0.099 


0.044 


2.28 


0.100 


0.056 


0.54 


0.030 


0.054 


0.62 


0.033 


1102.7 


3 


0.0144 3 4 


CR Boo 


0.027 


2.22 


0.060 


0.038 


1.59 


0.061 


0.058 


0.36 


0.021 


0.053 


0.45 


0.024 


1471.3 


5 


0.0109 6 


KL Dra 


0.049 


1.20 


0.059 


0.066 


0.91 


0.060 


0.094 


0.22 


0.020 


0.089 


0.27 


0.024 


1501.8 


7 


0.0193 7 


V803 Cen 


0.009 


5.32 


0.050 


0.017 


2.93 


0.051 


0.034 


0.54 


0.018 


0.030 


0.71 


0.021 


1612.0 


8 


0.00381 8 9 


CP Eri 


0.021 


2.17 


0.047 


0.035 


1.35 


0.048 


0.060 


0.28 


0.017 


0.054 


0.37 


0.020 


1701.2 


8 


0.00863 8 


V406 Hya 


0.017 


2.06 


0.034 


0.036 


0.99 


0.035 


0.076 


0.18 


0.014 


0.063 


0.26 


0.016 


2027.8 


10 


0.00673 11 



Skillman et~all jl999h 
Nelemans. Steeghs fc Grootl l|200lft 
Patterson et al.| |2002) 



O'Donoghue et al.l dl994h 
Provencal et al.l Jl997^ 



Patterson! lll997f) 
Wood et alJ j2002h 



Patterson! J200lf) 
Kato et al.l [|2004r i 



3 iRoelofs et all \200di ) 
k IWoudt fc Warned <2003l) 



to 



There als o exist s a form due to Eggleton that attempts to reconcile a low-mass appropriate form in Q3TJI with that due 
Nauenberg 1 19721 ) for higher masses 



0.0114 



/ M 2 \~s _ / Ah \ 



\M Ch J 
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3.5 



+ 



(33) 



where Mch = 1-44 is the Chandrasekhar mass and the constant M p = 0.00057 ( Verbunt fc Rappaport 19881 ) 



3 APPLICATION TO AM CVN SYSTEMS 



The results of applying the method outlined in section T2.2.2I to the AM CVn systems are shown in Table 1. Also given is 
the "dynamical only" result (ie. impress = 0) for comparison. The resu lting mass ratio for AM CVn is lower in every case 



than the normally quoted value q — 0.087 (eg. I Nelemans et al.l (|200lh ). as we have used the corrected, full expression for 
the dynamical expression (equation [4| . The effect of the addition of the pressure term is to require a larger mass-ratio to 
achieve the same precession rate. Since M 2 is largely determined by P or b, this in turn requires a smaller value for Mi. This 
effect notwithstanding, we see that the semi-degenerate TF model produces unacceptably high Mi values. In contrast, the 
degenerate relationship produces uncomfortably low values. 

Recent spectrosc opic observations o f the Hel 4471 line in AM CVn have yielded the first direct measurement of q for 
an AM CVn system I Roelofs et al. 20061 ). Disappointingly, their value of q — 0.18 is higher than any of those predicted by 
this superhump analysis. It is possible that the true effect of the pressure term for AM CVn is larger than that predicted 
from the above derivation but this begs the question of which unaccounted for factor has changed between the SU UMa and 
AM CVn systems? To achieve a result of q — 0.18 for the TF mass-radius relationship would require r/ne a factor ~ 7 larger 
than expected by the extension of the SU UMa expression. Examining the terms in equations (|27[) and (|30|l it is difficult to 



see any dependencies, even conspiring together, that could produce such a factor except cot i. The best fit value of rj H for 
SU UMa systems would correspond to a pitch angle i = 24° with the given fiducial parameters. This is sufficiently large that 
it may violate the tight- winding approximation that was used to derive equation ([3]). A systematic change in the value of 
i to 10° would produce the required factor. However, while this would give the "correct" result for AM CVn, a convergent 
solution that satisfies both the P OI h and P s h constraints simultaneously cannot then be found for any other system except 
HP Lib. The derived values for these two systems are listed in Table 2 calculated with the required scaling of uj pr for each 
mass-radius relationship to achieve an A M CVn mass ratio of q — 0.18. 

While simulations of low q binaries jTrussI 2007) do show a correlation between i and q, the trend is for i to increase 
as q decreases, ie. the exact opposite of what would be required to explain the change from SU UMa to AM CVn systems. 
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Table 2. Derived component masses after scaling the pitch angle i in the pressure term to give an AM CVn mass ratio 
q = 0.18 for each mass-radius relation. 



System 


TF (t = 10°) 
q Mi M 2 


SKH (i = 11°) 
q Mi M 2 


Degenerate (i 
q Mi 


= 15°) 
M 2 


Eggleton (i = 
q Mi 


15°) 
M 2 


AM CVn 
HP Lib 


0.180 0.649 0.117 
0.164 0.633 0.104 


0.180 0.517 0.093 
0.158 0.539 0.085 


0.180 0.186 
0.158 0.194 


0.033 
0.031 


0.180 0.206 
0.156 0.218 


0.037 
0.034 



Further, the same simulations show that the structure of accretion discs is largely unaffected by a significant decrease in the 
mass transfer rate that might be expected for systems with low-mass secondaries ruling out this as a possible cause. The 
simulations also show a change in the orientation of the accretion disc with respect to the line of centres of the two stars. This 
rotation angle 9 IC ,t appears to decrease with decreasing q for q < 0.1. The changes in i and rot may combine to explain why 
the discs in these simulations also show a tendency to be less centrally condensed at low q. In such a situation, the use of the 



Cannizzo &i Re iff (1992) profiles above may not be valid for such extreme systems. The effect appears to be sufficiently small 



however that it is not a viable explanation for the systematic change between the two classes and doubtful whether it is the 
explanation of the differences within the AM CVn group. 

A further factor that might be considered is whether the larger size of the primary's Roche lobe in extreme mass ratio 
systems might allow the disc to grow sufficiently that it could access the j = 2 resonance. Leaving aside the unique character 
of the j = 3 resonance that ena ble it to be ex cited l|Lubowlll991al ). evaluating u>dyn with j — 2 would actually cause the 
inferred q to be even smaller (cf. iPearsonl ([2003) that invoked resonances of higher j to explain a larger than expected mass 
ratio for AM CVn). Similarly, our analysis has characterised the disc response by the properties at a single radius. In reality, 
the precession is a collective property of the whole disc but perhaps we can find an effective radius r c a (not necessarily at a 
resonance) that gives good results? In summary, it is possible to achieve q = 0.18 for AM CVn with r e ff = and a degenerate 
mass-radius relationship. However, the value of Mi = 0.195 would then be unacceptably small. In contrast, r e ff = and 
the TF semi-degenerate secondary produces an acceptable Mi = 0.984 but a too small q = 0.118. We can satisfy both Mi 
and q simultaneously with r e g = rs and a TF semi-degenerate secondary that results in q = 0.19 and Mi = 0.607. However, 
applying this to other systems then produces unrealistically small values for Mi. These examples, th ough not exhaustive, 
strong ly suggest that it is not possible to find a suitable r e ff , M(R) combination. A study along the lines of lGoodchild fe Qgilvie 
1 2006) integrating across the whole disc structure would prove useful. 

The structure models are based on axisymm etric, vertically averaged but quite general accret ion disc equations with 
mean molecular weight and opacity as parameters ( Cannizzo fc Reifi 19921 ; Shakura fc Sunvaev 19731) . The critical transition 
temperatures are probably less well understood, although even here the arguments are sufficiently general that it is difficult 
to see that they could be so wrong as to explain the required change. We have tested the possibility that the proportionality 
coefficients were in error by trying to scale the results to give the correct result for AM CVn. Since this is not possible, the 
functional forms of the final equations themselves must be significantly in error. 

Probably the most crucial issue that should be considered is whether following the precedent of usin g a single mass-radius 



relatio n, such as those in equations (|31|l or (|32p . is appropriate for all the AM CVn secondaries. Figure 1 of lDelove. Bildsten fc Nelemans 



1 20051) 



shows that detailed models of the structure of the secondary produce a wide variation in radius for a given mass de- 
pending on the object's temperature. At the P or b of AM CVn, the majority of tracks are constrained to the region close to 
the fully degenerate relationship, but they increasingly diverge for smaller values of M%. The derived values of Mi in Table 2, 
however, seem more reasonable for the semi-degenerate forms of M<z{R%). In this case, the aforementioned figure would suggest 
a surprisingly hot secondary (logT 2 ~ 7.0-7.5). 

Regardless of the exact cause, the inability to find a universal coefficient for the strength of the pressure term that can 
be applied to all AM CVn systems forces us to conclude that there is no general relationship for these systems analogous to 
that found for the SU UMa binaries. The most likely reason is the lack of a single valued relationship between the secondary's 
mass and radius that makes it impossible to convert an expression involving P or b into one for M 2 . 



4 CONCLUSION 

We have shown how the extension of the theoretical relation relating P s h to q that works well for ordinary CVs appears to 
break down badly when applied to AM CVn systems. Hence, the values derived even by extrapolating the empirical relation 
from SU UMa to AM CVn binaries should be treated with extreme caution. The breakdown seems most likely to arise from 
the wide range of radii a secondary of a given mass can adopt. This may be exacerbated by a change in the accretion disc 
structure at very low mass ratios. Enhanced numerical studies that would enable such accretion disc profiles to be determined 
in terms of the system parameters may allow us to assess the importance of this contribution. Deriving a temperature for 
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such faint objects against the glare of the accretion disc would be very difficult, although in principle, it would enable us to 
place the secondary on an appropriate theoretical M2{Ri) track 
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